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Abstract 



We investigate the masses of + and 2 + glueballs in SU(2) lattice gauge theory using 
abelian projection to the maximum abelian gauge. We calculate glueball masses using 
both abelian links and monopole operators. Both methods reproduce the known full 
577(2) results quantitatively. Positivity problems present in the abelian projection are 
discussed. We study the dependence of the glueball masses on magnetic current loop size, 
and find that the + state requires a much greater range of sizes than does the 2 + state. 



1 Introduction 



The attempt to understand confinement in QCD has continued ever since Gell-Mann and 
Zweig introduced the quark concept. Attention has focussed almost universally on the 
pure gauge theory without light dynamical quarks as the place to begin study of the non- 
perturbative physics involved in confinement. In the absence of light dynamical quarks, 
what is meant by confinement can be precisely defined. Namely, in a confining theory, 
Wilson's loop must show an area law, or what is equivalent, the fundamental string tension 
must be finite. 

Having simplified full QCD to pure SU (3) gauge theory, many theorists have gone one 
step further and simplified to pure SU{2) gauge theory. The physics of confinement is 
thought to be quite similar for any of the SU(N) groups, of which the simplest is 577(2). 
While an SU(2) gauge theory obviously cannot explain the spectrum of baryons, lattice 
gauge theory simulations have shown that the meson spectra of pure SU(2) and 577(3) 
gauge theories are remarkably alike. 

A major effort has been devoted to the theoretical understanding of confinement in 
SU(2) pure gauge theory, both in the continuum and on the lattice. Within this frame- 
work, the dominant theme has been to search for a topological explanation of confinement. 
Topological objects are postulated to dominate the Euclidean path integral and lead to 
an area law for Wilson loops. The topological objects usually considered are instantons, 
monopoles, and vortices. Of these, only instantons are firmly established semi-classically 
in Euclidean pure SU(2) gauge theory. However, the general consensus of work done on 
instantons is that while they are surely highly relevant for understanding the breakdown 
of chiral symmetry in QCD, they appear to have little to do with confinement. 

This leaves monopoles and/or vortices as possible agents of confinement in pure SU(2) 
gauge theory. The present work is devoted to further exploring the monopole approach. 
Since the work of 't Hooft and Polyakov [|I|, monopoles are well-known to exist semiclas- 
sically if the gauge theory contains adjoint scalar fields in addition to the usual gauge 
fields. Such scalar fields are natural in N=2 supersymmetric SU(2) gauge theory, and 
the famous work of Witten and Seiberg showed that non-perturbative phenomena in this 
theory are indeed driven by the 't Hooft-Polyakov monopoles present in the theory |§. 
However, the relevance of the Witten-Seiberg work for understanding confinement in pure 
SU (2) gauge theory remains uncertain. 

The suggestion that monopoles nevertheless exist and control confinement in the pure 
non-abelian gauge theory was made much earlier by 't Hooft 0, who also advocated 
partial gauge-fixing as a way to see monopole physics more clearly. The basic physical 
picture of confinement via monopoles is the "dual superconductor", where the linear 
potential between heavy quark and anti-quark originates in a electric flux tube surrounded 
by the magnetic current of a monopole condensate. The fundamental string tension has 
been successfully calculated in this framework [|], The purpose of the present work 
is to see if the glueball spectrum can also be explained. Glueballs are as characteristic 
of a confining theory as the string tension, and any approach which purports to explain 
confinement should also produce quantitative results for the glueball spectrum. 
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To make the paper self-contained, the abelian projection formalism is reviewed in the 
next Section. The reader familiar with this can go immediately to Section We will 
work in lattice units throughout, and since all states are even under charge conjugation 
in S77(2) we will denote glueballs by their spin J and parity P as J p . 



2 The Abelian Projection 

In 't Hooft's framework, the first step is a partial gauge- fixing, applied only to those 
gauge fields which are "charged" , or have off-diagonal generators in the Lie algebra of the 
gauge group. The central idea is that the monopoles associated with the abelian gauge 
invariance left after partial gauge-fixing will control non-perturbative phenomena. For an 
SU (2) gauge group with the generator T 3 diagonal, the gauge field A 3 ^ is the abelian field 
or "photon", and gauge-fixing is done only on A a ,a = 1,2, or equivalently the charged 
fields 

After gauge-fixing, the charged field are ignored, the important physics being postulated 
to reside solely in the abelian gauge field and its monopoles. 

The 't Hooft approach has been explored rather extensively in SU{2) lattice gauge 
theory ||. The maximum abelian gauge has been found to be the only form of the gauge- 
fixing condition which allows a quantitative abelian calculation of the string tension. In 
the continuum form of the maximum abelian gauge, the continuum functional 

Gc^E/ (K) 2 + (4) 2 ) = Ww-)d*x (i) 

^ H J fi J 

is minimized over all SU(2) gauge transformations, leading to the conditions 

(d, + igA$W+ = (0„ - igAl)W- = 0, 

where g is the 577(2) gauge coupling. 

Stated precisely the purpose of the present paper is to see if the glueball spectrum can 
also be explained in this gauge. 



2.1 Lattice Gauge-Fixing 

The 577(2) lattice gauge theory is built out of link variables U^(x), 

where f = a /2 are the generators of SU(2) in the fundamental representation, and a is 
the lattice spacing. On the lattice, the maximum abelian gauge is obtained by maximizing 
the lattice functional 

tr r . 

(2) 
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over all SU(2) gauge transformations H. It is easy to show that in the continuum limit, 
maximizing Gi is equivalent to minimizing G c . 

The demand that Gi be stationary with respect to a gauge transformation at an 
arbitrary site y leads to the requirement that 

X{V) = E [UM^Uiiv) + Ul(y - fl)a 3 U(y - /})] (3) 

be diagonal. This can be accomplished by a gauge transformation However, the 

value of X at the nearest neighbors of y is affected by fl(y), so the diagonalization of X 
over the whole lattice must be done iteratively. 

It is useful to factor a U{1) link variable from U^(x), writing U^(x) = u^(x)w^(x), 
where u^(x) = exp(z<^r 3 ), and w^{x) = exp(^ ■ r), with 6^ = 0. The angle 0^ can be 
extracted from the matrix elements of by expanding the gauge-fixed SU(2) link in 
Pauli matrices, writing 

U» = U° + iJ2K-cr k . (4) 

k=l 

Then 0^ = 2 • arctan(£/^/£/°). The SU(2) action is periodic in 0^ with period An. To 
carry over standard formulas from U(l) lattice gauge theory, it is more convenient to use 
an angle 0j = 0J/2, with 0^ G (-tt, vr]. 



2.2 Monopole Location in SU(2) 

The location of the magnetic current starts with plaquette angles (jp constructed from 

$U*) = 9^1 - d v 4>l = 4>l{x) + 4>l{x + //) - 0j(x + u)- 4>l(x) (5) 

The plaquette angle (j? is resolved into a Dirac string contribution, plus a fluctuating 
part: 

4>%, = 27r *V + 4>%,, (6) 

where ffi G (— 7r,7r] and is an integer ||. The integer- valued magnetic current m M is 
determined by the net flux of Dirac strings into an elementary cube: 

= --e^p ■ d v n a p. (7) 

In the identification of m^, elementary cubes were used, so a magnetic charge is located 
at the center of a spacial 1-cube, likewise for other components of the magnetic current. 
This is done mainly for practical reasons; if (say) 2-cubes were used instead, the effective 
lattice size would become 8 4 instead of 16 4 . 
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2.3 Abelian Wilson Loops 

The maximum abelian gauge has in a global way made the charged link angles 9^ as 
small as possible. Non-perturbative quantities are postulated to be solely contained in 
the abelian link angles 0^ or equivalently (fp . This allows two abelian approximations 
to the full SU(2) Wilson loop. In the first, we simply use the link angles in the 

standard U(l) formula for a Wilson loop: 

Wu(i) = (exp (i 4>lJ» 

where (■) denotes the expectation value over the ensemble of 3 link angle configurations, 
and Jfj, is an integer- valued current defined by the path of the heavy quark. 

In £7(1) lattice gauge theory, a Wilson loop expressed as in Eq.(§) can be factored into 
a term arising from photon exchange, times a term arising from monopoles, 

W U{1) = W phot ■ W mon . (9) 

The photon contribution is perturbative, so for non-perturbative quantities, the same 
information should be in W mon as Wua)- The monopole Wilson loop W mon forms our 
second abelian approximation to the full £77(2) Wilson loop. An explicit formula for 
W mon is obtained by writing J M as the curl of a Dirac sheet variable 0; = d^D^, 
where d v denotes a discrete derivative. Then W mon is given by 



W mon =(exp — Y^D^FUx) )) , (10) 




where {■) m denotes the sum over configurations of magnetic current. For a rectangular 
Wilson loop, a useful choice for the sheet variable D^ v is to set D^ v = 1 on the plaquettes of 
the flat rectangle with boundary J M , and D^ u = on all other plaquettes. In Eq. fllOD , F* 
is the dual of the field strength due to the magnetic current; F*(x) = ^€^ ua pF a p(x). The 
field strength itself is derived from a magnetic vector potential A™, F^ v = d^A™ — d u A™, 
where 

A™(x) = Y.v(x-y)m»(y), (11) 



y 



and fh/j, is the integer- valued, conserved magnetic current defined in Eq.(|7|). 

To summarize, once we have gone to the maximum abelian gauge, our two abelian 
approximations to SU{2) Wilson loops are Eq.@ or Eq. (|T0D . If the basic assumptions 
involved in the abelian projection are correct, non-perturbative quantities like the string 
tension and glueball spectrum should be obtainable from either one. 



lr rhere is in principle a finite volume correction to Eq. (|To|) arising from Dirac sheets, fiofl , but we have 
checked that it is negligible in the present calculations. 
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3 Glueball operators 



We will make use of standard methods for extracting the glueball spectrum in lattice gauge 
theory, as expounded in particular by Teper and collaborators |Tl|]. The basic object of 
interest is $, called the glueball "wave function". This is a gauge-invariant combination 
of SU(2) links constructed so as to make as large as possible, where \Q > is 

the SU(2) vacuum, and \G > is a one glueball state, usually taken to have a definite 
three-momentum. In practice, $ is a relatively simple Wilson loop-like combination of 
SU(2) links, summed over orientations to project a definite spin-parity, and over spacial 
locations to project a definite three-momentum. An important advance made by Teper 
was to notice that | < > | is numerically much larger if $ is constructed out of 

"fuzzy" links [12||. Fuzzy links are constructed by replacing original products of two links 
in a given direction by themselves plus some fraction of "staples" (four-link products with 
the same endpoints as the original link). A related procedure applied directly to single 



links was introduced by the Ape collaboration [13|. This is commonly called "smearing' 



and is what we use in the present work, when calculating with abelian links. To keep 
everything strictly local in Euclidean time, the staples used in smearing are spacial, i.e. 
for a z— link, the smeared staples are x — z — x and y — z — y, but not t — z — t. Smearing 
is typically done several times in the full gauge theory ||13|| , but here we smear the abelian 
links once, twice, or three times. 

For the + and 2 + glueballs, $ is constructed out of simple rectangular Wilson loops, 
anchored at the spacial point x, and lying in a particular spacial plane. The spin projection 
for + is done by averaging over the three spacial planes, x — y, x — z, and y — z, while for 
the 2 + , the spin projection is accomplished by subtracting twice the x — y rectangular 
loop from the sum of z — x and z — y loops ||11|| . Projection to a three- momentum p = 
is done by summing over spacial points x, divided by N x N y N z , where N x , N y , and N z are 
the number of lattice sites in x,y, and z-directions. 



3.1 Glueballs and the Abelian Projection 

In the present work, we evaluate the Wilson loops in the wave functions $ either in terms 
of the abelian link angles 3 as in Eq.(^|), or in terms of the magnetic current of monopoles 
as in Eq . (|T0|) . Using these, we calculate correlation functions which we label as 

C(t) u{1) = ($(t)$(0)) l/(1) /($(0)$(0)) f/(1) 

C(t) mon = ($(t)$(0)) mo „/($(0)$(0)) mon , 

where t refers to the separation between the two time-slices where $ is evaluated. Any 
vacuum expected value of $ is assumed to have been removed. These calculations are 
analogous to calculating the potential using Eq.(||) and Eq.(|10]), respectively. The aim 
of course is to see if the glueball spectrum can be obtained in the maximum abelian 
gauge to a similar degree of accuracy as the string tension. The use of smeared links 
is straightforward in the calculation which uses abelian links directly, namely C(t)j/(i), 
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and for this case we will report on our results both with and without smeared abelian 
links. The Wilson loops in the monopole correlation function are obtained by computing 
the dual flux through the surface of the loop. Calculating this flux through loops with 
smeared edges is somewhat complicated, so for the monopole correlation function we 
simply calculate directly using the simplest plane rectangular surface . In any case, after 
gauge-fixing to the maximum abelian gauge, smearing is not essential for obtaining good 
results, unlike the situation for full SU(2). The process of going to the maximum abelian 
gauge and using either the resulting abelian links or monopoles results in correlation 
functions with much less noise than the full 577(2) correlation functions. 

3.2 Positivity 

The full SU{2) correlation function C(t) S u(2) is of the form 

C{t) sm = S,c 4 ■ e~ Eit (12) 

where we again assume any vacuum expected value of $ has been removed, and for 
simplicity treat the case of a system finite in space but infinite in time. Since C(t)su(2) is 
normalized to 1.0 at t = 0, we have 

Q = iw>r/£iw>r> 

i 

so the non-zero coefficients q are positive and certainly less than unity. On a lattice 
with finite time extent and periodic boundary conditions, we must make the replacement 
e -Eit _^ e -Eit _j_ e -Ei(T-t) j n Eg.( |l2]) to satisfy periodicity. In addition, a constant b, 
independent of t but of order exp —T will in general be present in the correlation function 
|14j| . It is easy to show that for an SU(2) action with nearest neighbor couplings in time 



and a positive transfer matrix, the coefficients c% continue to be positive and < 1.0. Our 
calculations were done with the Wilson 577(2) action, which satisfies the above properties. 

The correlation function C(t)um is obtained by replacing each 577(2) link in $ by its 
£7(1) approximation in the maximum abelian gauge. The correlation function C(t) mon in- 
volves a further approximation in which only the contributions of monopoles are retained, 
and photons or abelian gluons are dropped. 



The correlation functions C(t)u(i) and C(t) mon are fitted to a form [14 



6 + c -(e-^ + e-^ T -*)). (13) 

The fit is aimed at getting only the lowest energy eigenvalue, Eq, in a given channel. The 
results show that the constant b is always very small, essentially zero, while the value of 
of Eq is generally in agreement within statistical errors with the corresponding full SU(2) 
value. 

Problems with positivity show up in two ways. The first has to do with the constant 
Co- In the full 577(2) correlation function a good choice of the operator $ will generally 
increase the value of Co, but positivity always requires Co < 1.0. However in C{t) mon and 
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C(t)u(i), fits of the form of Eq.(^) do result in Co > 1.0. The problem is more severe for 
C(t) mon where values of c > 1.3 are common. 

The second difficulty is that an acceptable fit of the form Eq.fllBD is often not possible 
without dropping small t points, t = 0, t = 0, 1, or in some extreme cases t — 0,1, 2. This 
could happen for a function which does satisfy positivity, if the lowest eigenvalue E Q is 
not sufficiently dominant. In such a case, the 'effective eigenvalues' 

E {t) = -ln{C{t + l)/C{t)) 

will decrease with increasing t, reflecting the fact that the exact correlation function of 
the form Eq.(^), is concave upward. The failure of the fit to a simple exponential would 
in this case merely represent the importance of physical excited states in the correlation 
function. However for C(t)u(i) and C(t) mon , it is often true that the values of E (t) 
increase with t for t — 0,1, showing that C(t)u(i) and C(t) mon violate this concavity 
requirement for small values of t. Our strategy for dealing with this problem violation is 
simple. We examine the effective eigenvalues E (t) and look for a 'plateau'. If for small 
values of t, the E (t) are less than those in the plateau, those values of t are omitted in 
the fit. The fit then runs from some t m i n to t=7 on our 16 4 lattice. This procedure is 
of course heuristic, and is based on the assumption that for sufficiently large values of t, 
the true lowest eigenvalue of the underlying SU(2) gauge theory will dominate in C{t)u{i) 
and C(t) mon . 

In a given channel, we use rectangular Wilson loops of various sizes as operators, and 
compute only diagonal terms, i.e. correlations of the same operator with itself. In the 
resulting fits, we favor results where the value of cq is as close to 1.0 as possible, that is 
we disfavor small Co since it implies that the operator has a small overlap with the state 
of interest, and we also disfavor values of Cq much larger than 1.0 since they imply a large 
amount of positivity violation. Detailed examples will be given in the next section, where 
it will be seen that despite these problems with positivity, the final results are surprisingly 
reliable. 

It is perhaps worth making some general comments as to why positivity violation 
occurs at all. Since we have a gauge system, a configuration of link angles on a given 
time-slice represents a physical configuration plus many redundant variables. One of the 
reasons for gauge-fixing is to reduce the number of irrelevant variables. However, as 
opposed to a gauge which acts on some local operator, putting a link angle configuration 
into the maximum abelian gauge is an iterative process, involving many sweeps of the 
lattice, and therefore coupling sites which are far apart on the lattice. While the non- 
locality thereby introduced must cancel out if only SU (2) gauge-invariant operators are 
used, this is not true when links are represented only by their U(l) parts. The effective 
action which generates U(l) configurations in the maximum abelian gauge no longer has 
the properties of the original Wilson SU(2) action, in particular, it will not be nearest- 
neighbor in time. Correlation functions can then violate standard positivity properties. 
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4 Calculations 



4.1 Simulation and Gauge-Fixing 

Our simulations were done on a 16 4 lattice, using the standard Wilson form of the SU(2) 
action, at (5 — 2.40, and 2.50. At each /3, after equilibration, 2000 configurations were 
saved. Saved configurations were separated by 20 updates of the lattice, where a lattice 
update consisted of one heatbath sweep fL5f , plus one or two overrelaxation sweeps [|16j. 



Each of these configurations was then projected into the maximum abelian gauge using 
the overrelaxation method of Mandula and Ogilvie, with their parameter uj = 1.70 ||17|| . 
The overrelaxation process was stopped after the off-diagonal elements of X(x) of Eq.((^) 
were sufficiently small. Expanding X(x) in Pauli matrices, 

3 

X = X\x)+iY,X\x)-a kl (14) 
fc=i 

we used 

(|x c Y) = ^E(l*W + |x 2 (*)| 2 ) (is) 

as a measure of the average size of the off-diagonal matrix elements of X over the lattice, 
and required (jX ch \ 2 ^ < 10~ 10 . This condition was reached in approximately 1000 over- 
relaxation sweeps. Even though our condition on ^|X c/l | 2 ^} is very tight, the gauge-fixed 
configuration is not unique; Gribov copies can still occur. However, Hart and Teper have 
shown that the uncertainty in the string tension caused by the Gribov ambiguity is rather 



small ||18|| . We have assumed the same will hold for the glueball spectrum. 

From each gauge-fixed SU(2) link, the U(l) link angle 0^ was extracted using the 
formula = arctan([/^/£/°), as described in Section [[]. Then making use of the plaquette 
angles ({P , the magnetic current fh^{x) was found. 



4.2 The Glueball Spectrum in Maximum Abelian Gauge 

In this section, we present our results for the glueball spectrum in detail. Our fits are 
all of the form Eq. (|T3D . In Tables 1-8 below we present the data from our fits for the + 
and 2 + channels at zero momentum for for abelian links and monopoles at (3 = 2.50, and 
(3 = 2.40. Since the three-momentum is zero, Eq is in fact an estimate of the mass of 
the corresponding glueball , denoted as M in the tables. The column labelled 'operator' 
labels the Wilson loop used in calculating C(t)u^ or C{t) mon . Thus in Table 1, the 
operator 2x1 means the correlation function involves a 2 x 1 Wilson loop correlated with 
itself, while the operator (2 x means a 2 x 1 loop smeared once, etc. The column 
labelled 'fit' lists the small t points omitted, if any, in the fit. The next column gives the 
fitted value of the coefficient Co in the fit, which we henceforth call the 'overlap'. 

We have also studied (less extensively) the glueballs masses as derived from correlation 
functions with finite momentum. We took the momentum to be the lowest possible on 



9 



the lattice, 7r/8 in lattice units, in (say) the ^-direction. The value of E was derived from 
the corresponding correlation function in the same manner as for zero momentum. Then 
the mass M was obtained by assuming the usual (continuum) relation between E and 
M , 



E = VM^+p 2 . 

Representative finite momentum results at (3 = 2.50 are presented in Table 9. 

Before presenting our final estimates of glueball masses, we first remark on some 
systematic trends which can be seen in the various tables. For data obtained from the 
link correlation functions (Tables 1-4), the overlap coefficient Cq is safely within error bars 
of the physical range < Cq < 1 for almost all cases. In contrast, for the monopole 
correlation functions (Tables 5-8), Cq > 1 is the rule rather than the exception. Here we 
only find cq < 1 for the first few operators for the 2 + state at (3 = 2.50. So for the link 
correlation functions, this gross violation of positivity is rather uncommon, while it is 
quite common for the data obtained from monopole correlation functions. Nevertheless, 
the values of Eq given in the tables for the two cases indicate that link and monopole 
operators are coupling to the same set of states. 

The need to drop small values of t from the fits to Eq.(^) is another indication of 
positivity violation, caused in each case by small values of the effective eigenvalues E (t). 
This dropping of some points is always necessary when cq > 1, the extreme case being 
the 2 + at (3 = 2.5 state from monopoles, where c > 1.4 required dropping t = 0,1,2. 
However, for c < 1, low effective eigenvalues for small t also necessitate dropping points 
at small t. The 2 + at (3 = 2.5 state from monopoles also illustrates this. 

We now turn in Table ^ to a comparison of our results for glueball masses with the 
best available full SU (2) data. The latter is from Michael and Teper |19[ who extracted 



glueball masses from 20,000 measurements at /3 = 2.40 on a 16 4 lattice, and 14,000 at 
(3 = 2.50 on a 20 4 lattice. Generally the agreement is good between monopoles and links 
and with the full SU(2) results, especially for the lightest and heaviest states, namely 
the + at (3 = 2.50 and the 2 + at (3 = 2.40. For the remaining states, there are some 
discrepancies at roughly the 5% level. This is not unreasonable, given the heaviness of 
glueballs vs. the square root of the string tension, and the relatively small number of 
measurements (2000) used in the present work. 



4.3 Magnetic Current Loop Size and the Glueball Spectrum 

In the maximum abelian gauge, the magnetic current can be resolved into closed loops 
containing different numbers of links. At (3 = 2.50 , the percentage of magnetic current in 
loops containing less than 10, 20, 50 and 100 links is 43%, 51%, 56%, and 59%, respectively 
[|J. It has been known for some time that the fundamental string tension is unaffected 
by the small loops of magnetic current. At (3 = 2.50, dropping all current contained in 
loops with less than 100 links does not affect the string tension. On the other hand, if we 
increase the size of the smallest allowed magnetic current loop to 200 links, the resulting 
string tension is approximately 10% smaller than the correct value. So we can say that 
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loops ~ 100 links and larger are necessary to explain the (3 = 2.50 string tension. 

In this section, we explore the dependence of the glueball spectrum on magnetic current 
loop size, discussing for simplicity only f3 = 2.50. For the + state, we will use 2x1 Wilson 
loops in C(t) mon , for the 2 + state, we will use 2x2 Wilson loops. The first question we 
asked is whether the masses are stable under the elimination of small magnetic current 
loops. We tested this by successively eliminating all loops with less than 10, 20, 50, 100 
etc. links. If the glueball masses were to behave in the same manner as the string tension, 
the masses would be preserved up to a lower limit on sizes of ~ 100 links, but then would 
fall below their correct values. In fact the glueball masses behave somewhat differently 
than the string tension. For the + state, using the entire magnetic current, and the 
2x1 Wilson loop in C(t) mon , we obtained a mass of 0.72(2). While we obtain values 
consistent with this when we eliminate all magnetic current loops smaller than 6, 10, and 
20 links, by the time we eliminate all loops containing less than 50 links , the + mass 
has fallen to 0.60(3) , and steadily declines as we increase the size of the smallest loop 
of magnetic current allowed. Thus the + glueball state is sensitive to smaller loops of 
magnetic current than the string tension. 

For the 2 + state, using the entire magnetic current and the 2x2 Wilson loop in 
C(t) mon , we obtained a mass of 0.98(8). We find that this mass is stable in successively 
eliminating loops of size less than 10, 20, 50, 100 links , and even somewhat beyond. (For 
the lower limit of 100 links, we obtained a mass of 0.98(2).) So the lower bound on sizes 
of loops of magnetic current needed to reproduce the 2 + glueball is quite similar to that 
needed for the string tension. 

Since glueball masses, like the string tension, are non-perturbative quantities, it might 
seem obvious that they would require the largest loops of magnetic current. The presence 
of these very large loops of magnetic current is thought to be a defining characteristic 
of a confining phase. For example, in U(l) lattice gauge theory, they are present in 
the confining phase, and absent in the deconfined Coulombic phase p0| |. To investigate 
this in the SU(2) glueball spectrum, we tried extracting the masses from C(t) mon after 
eliminating all loops of magnetic current larger than a certain size. We find that the + 
and 2 + states behave rather differently under these cuts. For the + , if the upper limit on 
loop sizes is taken to be 100 links, the resulting mass is close to the value obtained with 
the entire magnetic current, 0.75(3), vs. 0.72(2). However, upon raising the upper limit 
to 200 links, the mass is too small , 0.53(5) , and remains so when the upper limit on loop 
size is raised still further giving 0.63(5) when the upper limit is 1000 links. Only when 
the upper limit is raised greater than 1000 links does the mass finally stabilize, reaching 
0.70(8) by the time the largest loop of magnetic current contains 1500 links. Thus the + 
state is sensitive to the largest loops of magnetic current. 

For the 2 + , the story is different. When the upper limit on loop size is 100 links, the 
mass value obtained is too large, 1.48(6), vs. 0.98(8) obtained with the entire magnetic 
current. But by the time the upper limit is raised to 200 links, the mass is 1.05(15), and 
remains within error bars of 0.98(8) as the upper limit is raised steadily to 2000 links. So 
unlike the + and the string tension, explaining the 2 + mass does not appear to require 
the largest loops of magnetic current. 
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In Table [II], we show + and 2 + mass at f3 — 2.50 with various cuts on the number 
of links allowed in loops of magnetic current. In the Table, N min denotes the number of 
links in the smallest allowed magnetic current loop, and N max does likewise for the largest. 
There are two rather intriguing results. The first is the + state basically requires all of 
the magnetic current which resides in loops > 50 links in size. (Note that this still means 
that over 50% of the current can be discarded.) The second is the relative insensitivity of 
the 2 + state to the particular cut made on magnetic current. The 2 + state is visible in a 
wide variety of windows, the exception being one which only allows small loops. 

An interesting difference between + and 2 + states has also been found in the work 
of Schafer and Shurak [[H]] on instantons and glueballs. We now discuss the possible re- 
lation of their results to our own. In the instanton case, the + channel receives strong 
non-perturbative contributions from individual instantons and anti-instantons, as well as 
instanton-anti-instanton interactions. The 2 + channel, on the other hand, does not re- 
ceive contributions from individual instantons and anti-instantons, only a rather weak 
contribution from instanton-anti-instanton interactions. On the lattice, it is known that a 
distribution of instantons and anti-instantons, when cast into the maximum abelian gauge, 
is represented by a network of magnetic current [22], and a rough correspondence can be 
made between the size of an instanton and the size of the magnetic current loop it gener- 
ates. Since instantons, anti-instantons and their interactions all give a non-perturbative 
contribution to the + correlation function, but only the interactions contribute to the 
2 + , it is natural for the + state to depend more strongly on magnetic current loops of all 
sizes than the 2 + . While this conclusion is similar in the two approaches, our work would 
imply that the + still could not be quantitatively explained in an instanton gas/liquid, 
since the latter, having no confinement |23|, contains too few of the very largest loops of 
magnetic current. 



5 Conclusions and Summary 

We have shown that contact can be made with the spectrum of glueballs in SU(2) lat- 
tice gauge theory, using the maximum abelian gauge and computing correlation functions 
using either abelian links or monopoles. The calculations we have presented show rather 
convincingly that for sufficiently large values of t, both C(t)u^ and C(t) mon are domi- 
nated by the lightest glueball in the channel of interest. The present calculations were 
done starting with 2,000 SU(2) configurations, as opposed to the O(20, 000) used in con- 
ventional SU{2) calculations. In addition, when glueballs are found using C mon (t), only 
the magnetic current is needed, which occupies around 1% (at (3 = 2.50) of the links on 
the lattice [|[]. So despite the effort required to project a configuration into the maxi- 
mum abelian gauge, the description finally involves a rather small subset of the variables 
in a full SU(2) configuration. Further, the results on how different components of the 
magnetic current build up the + and 2 + give information not readily obtainable in a 
conventional SU(2) calculation. 

Nevertheless, the present method of calculation; gauge-fixing, followed by a trunca- 
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tion of the degrees of freedom, is not a controlled approximation, and the results must be 
checked by comparing them with those from full £77(2) calculations. One of the major 
uncertainties which shows up in the glueball calculations reported here is the presence of 
positivity violation. In calculations of the heavy quark potential, this was hardly visible 
at all. The only sign of it was a very small wrong sign Coulomb term seen in the potential 
calculated from monopoles [|J. No positivity violation was seen in the potential calcu- 
lated with abelian links. In contrast, the glueball calculations done here show positivity 
violation occasionally with abelian links, and quite commonly with monopoles. It would 
clearly be of great interest to know the spectrum of masses which are allowed to appear 
in C(t)u(i) and C(t) mon . If only physical states can occur, then positivity violation would 
imply that excited states may have coefficients with negative signs, i In this situation, 
standard upper bound statements on glueball masses are lost, but the spectrum still con- 
tains only physical states. In the maximum abelian gauge, the greatly reduced statistical 
errors of C(t)u(i) and C(t) mon vs C{t) S u(2) might make this a price worth paying. 

We are interested in extending the present work in several directions. The first goal is 
to include the 0~ and 2~, the next states up in mass in the SU(2) glueball spectrum. Be- 
yond this, we want to gather a much larger dataset to pin down any possible discrepancies 
between abelian link/monopole glueball masses and those obtained with full SU(2), and 
finally attempt to extract excited state masses and compare those with full SU (2) . This 
could shed light on the allowed spectrum of C(t)u(i) and C(t) mon , discussed above. The 
explanation of the glueball spectrum is clearly a challenge that any topological approach 
to confinement must meet. 

This work was supported in part by the National Science Foundation under Grant 
No. NSF PHY 94-12556. During the early stages of this work, R. F. was supported in 
part by Department of Energy grant DE-FG02-91ER40677. The calculations were carried 
out on the Cray T-90 system at the San Diego Supercomputer Center at the University 
of California, San Diego, and the Power Challenge system at the National Center for 
Supercomputing Applications at University of Illinois, both supported in part by the 
National Science Foundation. 
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Tables 



Table 1: + state from abelian links at (5 — 2.50 



operator 


fit 


Co 


M 


2 x 1 





0.81(2) 


0.76(2) 


(2 x 1)_! 




0.99(1) 


0.68(1) 


2x2 





0.91(2) 


0.70(2) 


(2 x 2)_x 


0,1 


1.15(8) 


0.68(4) 


3x2 





0.92(2) 


0.67(2) 


3x3 


0,1 


1.04(7) 


0.67(4) 


4x4 


0,1 


1.04(7) 


0.63(4) 



Table 2: 2 + state from abelian links at (3 — 2.50 



operator 


fit 


Co 


M 


3x3 





0.63(2) 


1.17(5) 


(3 x 3)_x 





0.95(3) 


1.06(3) 


(3 x 3)_ 2 




0.992(5) 


1.04(1) 


4x4 





0.62(6) 


1.02(4) 


(4 x 4)_! 





1.05(3) 


0.99(3) 


5x5 





0.55(2) 


0.94(4) 



Table 3: + state from abelian links at (5 — 2.40 



operator 


fit 


c 


M 


1 x 1 





0.93(3) 


1.10(4) 


(1 X 1)_! 




0.94(3) 


0.99(3) 


2 x 1 





0.99(3) 


1.07(3) 


2x2 


0,1 


1.2(2) 


1.04(7) 
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Table 4: 2+ state from abelian links at (3 — 2.40 



operator 


fit 


c 


M 


2x2 





0.67(9) 


1.7(1) 


(2 x 2)_i 





0.86(8) 


1.5(1) 


3x3 


0,1 


0.73(6) 


1.48(9) 


(3 X 1)_! 





0.7(1) 


1.6(1) 


(4 x 2)_! 




1.00(1) 


1.45(2) 



Table 5: + state from monopoles at (5 — 2.50 



operator 


fit 


c 


M 


2x1 





1.19(2) 


0.72(2) 


2x2 


0,1 


1.38(9) 


0.70(4) 


3 x 1 





1.19(2) 


0.69(2) 


3x3 


0,1,2 


1.9(4) 


0.72(8) 


4x4 


0,1,2 


2.0(4) 


0.67(7) 



Table 6: 2 + state from monopoles at (5 — 2.50 



operator 


fit 


c 


M 


2x2 


0,1 


0.9(1) 


0.98(8) 


3 x 1 


0,1 


0.6(1) 


1.0(1) 


4 x 1 


0,1 


0.7(1) 


0.93(8) 


5 x 1 


0,1 


0.7(1) 


0.94(9) 


3x2 





1.21(3) 


0.95(2) 


3x3 


0,1,2 


1.4(5) 


0.9(1) 


4x4 


0,1,2 


1.8(6) 


0.9(1) 



Table 7: + state from monopoles at (3 — 2.40 



operator 


fit 


c 


M 


1 x 1 





1.33(4) 


1.04(2) 


2 x 1 





1.34(3) 


0.98(2) 


2x2 


0,1 


1.6(1) 


0.97(5) 


3 x 1 





1.34(3) 


0.94(2) 
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Table 8: 2 + state from monopoles at (5 — 2.40 



operator 


fit 


Co 


M 


2x2 





1.39(7) 


1.39(6) 


3 x 1 





1.35(7) 


1.53(5) 


4 x 1 





1.25(6) 


1.47(5) 


5 x 1 





1.23(8) 


1.48(7) 


6x1 





1.17(8) 


1.47(7) 



Table 9: Finite momentum results from links and monopoles at (3 — 2.50 



state 


type 


operator 


Co 


M 


0+ 


links 


2x1 


0.81(2) 


0.74(3) 


2+ 


links 


4x4 


0.66(3) 


1.06(5) 


0+ 


monopoles 


2 x 1 


1.21(2) 


0.66(2) 


2+ 


monopoles 


2x2 


1.24(4) 


1.07(3) 



Table 10: Glueball masses from full SU(2), links, and monopoles 



p 


state 


SU(2) 


links 


monopoles 


2.5 


0+ 


0.72(3) 


0.69(1) 


0.70(1) 


2.5 


2+ 


1.05(3) 


1.04(1) 


0.95(2) 


2.4 


0+ 


0.94(3) 


1.06(2) 


0.99(1) 


2.4 


2+ 


1.52(3) 


1.47(2) 


1.49(4) 



Table 11: The + and 2 + masses with cuts on the magnetic current 



N • 


N 

1 ' max 


0+ 


2+ 





OO 


0.72(2) 


0.98(8) 


20 


OO 


0.69(4) 


1.00(2) 


50 


OO 


0.61(3) 


0.97(2) 


100 


OO 


0.49(2) 


0.98(2) 


20 


2000 


0.69(5) 


0.98(2) 


20 


1000 


0.51(5) 


1.00(3) 


50 


2000 


0.61(5) 


0.98(2) 


100 


1000 


0.40(6) 


0.96(3) 


100 


500 


0.1(1) 


0.95(4) 
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Figure Captions 



Figure 1. The + correlation function from abelian links. Symbols: 
diamonds — (2x1)! loops, crosses — 2x2 loops , squares —3x2 loops. 

Figure 2. The + correlation function from monopoles Symbols: 
diamonds — (2 x 1) 1 loops, crosses —2x2 loops , squares —3x1 loops. 

Figure 3. The 2 + correlation function from abelian links. Symbols: 
diamonds — (3 x 3)i loops, crosses — (4 x 4) 1 loops , squares —4x4 loops. 

Figure 4. The 2 + correlation function from monopoles. Symbols: 
diamonds —2x2 loops, crosses —3x2 loops , squares —3x1 loops. 
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